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Spectral energy distribution of the 7-ray microquasar LS 5039 
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Abstract. The microquasar LS 5039 has recently been detected as a source of very high energy (VHE) 7-rays. This 
detection, that confirms the previously proposed association of LS 5039 with the EGRET source 3EG J1824— 1514, 
makes of LS 5039 a special system with observational data covering nearly all the electromagnetic spectrum. In 
order to reproduce the observed spectrum of LS 5039, from radio to VHE 7-rays, we have applied a cold matter 
dominated jet model that takes into account accretion variability, the jet magnetic field, particle acceleration, 
adiabatic and radiative losses, microscopic energy conservation in the jet, and pair creation and absorption due to 
the external photon fields, as well as the emission from the first generation of secondaries. The radiative processes 
taken into account are synchrotron, relativistic Bremsstrahlung and inverse Compton (IC). The model is based 
on a scenario that has been characterized with recent observational results, concerning the orbital parameters, 
the orbital variability at X-rays and the nature of the compact object. The computed spectral energy distribution 
(SED) shows a good agreement with the available observational data. 
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1. Introduction 

Microquasars are radio emitting X-ray binaries with rel- 
ativistic radio jets (Mirabel & Rodriguez [1555) 1 . LS 5039 
was classified as microquasar when non-thermal radiation 
produced in a jet was detected with the VLBA (Paredes 
et al. I2000[) . New observations conducted later with the 
EVN and MERLIN confirmed the presence of an asym- 
metric two-sided persistent jet reaching up to ^1000 AU 
on the longest jet arm (Paredes et al. [2002 . Rib6 l2552|l . 
LS 5039 is a high mass X-ray binary (Motch et al. I1997f) 
whose optical counterpart is a bright (V <~ 11) star of 
spectral type Q 6.5V((f)) (C lark et al. 125511) . Although 
McSwain et al. I|2001[ 12554) first derived the period and 
the orbital parameters of the system, new values of the or- 
bital period (P or b = 3.9060 ± 0.0002 days, with periastron 
passage associated to T =HJD 2 451 943.09±0.10), eccen- 
tricity (e = 0.35 ±0.04), phase of periastron passage (0.0), 
and inclination (assuming pseudo-synchronization, i = 
24.9 ± 2.8°) have been reported by Casares et al. 
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who place the system at a distance of 2.5 ±0.1 kpc. The 
radius of the companion star is moderately close to the 
Lagrange point during periastron but without overflowing 
the Roche lobe. In addition, Casares et al. 1|2005[) state 
that the compact object in LS 5039 is a black hole of 
mass Mx = 3.7l};o M . At X-rays, the spectrum of the 
source follows a power-law and the fluxes are moderate 
and variable, with temporal scales similar to the orbital 
period. The variability is likely associated with changes in 
the accretion rate due to the motion of the compact object 
along an eccentric orbit (Bosch- Ramon et al. 12505a). 

Among microquasars, LS 5039 is specially relevant for 
its proposed association with 3EG J1824— 1514 (Paredes 
et al.[2p00), an unidentified source in the 3rd EGRET cat- 
alog of high-energy 7-ray sources (>100 MeV; Hartman 
et al. H999|) . Recently, Aharonian et al. (2005a), using the 
High Energy Stereoscopic System (HESS), have detected 
LS 5039 at energies above 250 GeV. This detection, that 
confirms the high-energy 7-ray nature of LS 5039 (Paredes 
et al. I2000|) . gives strong support to the idea that mi- 
croquasars are a distinctive class of high- and very high- 
energy 7-ray sources. The confirmation of LS 5039 as a 
high-energy 7-ray source, combined with the fact of being 
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a runaway microquasar ejected from the Galactic plane at 
~ 150 km s _1 (Ribo et al. 12002(1 . strongly raises the ques- 
tion about the connection between some of the remain- 
ing unidentified, faint, variable, and soft 7-ray EGRET 
sources and possible microquasar systems above/below 
the Galactic plane (e.g., Kaufman Bernado et al. [2002 
Romero et al.HUl Bosch-Ramon et al. l2005bjl . 

LS 5039 was also proposed to be the counterpart of a 
COMPTEL source (-1 MeV) (Strong et al. l2TO]l and 
BATSE detected it at soft 7-rays (Harmon et al. I2004J1 . 
In addition to LS 5039, there is another microquasar, 
LS I +61 303 (Massi et al. I2UU4) . that is likely associ- 
ated with an EGRET source (Kniffen et al. H557) l. Both 
microquasars have been extensively studied at different 
wavelengths, and both have a very similar spectral energy 
distribution. Bosch-Ramon & Paredes (2004a, b) have ex- 
plored with a detailed numerical leptonic model whether 
these systems can produce the level of emission detected 
by EGRET (MOO MeV), an d the o bserved variability. 
At VHE, Aharonian et al. {2005b) argue in favor of 
hadronic origin of TeV photons in case that 7-rays are pro- 
duced within ~ 10 12 cm, and Dermer & Bottcher (2005) 
point to combined star IC and synchrotron self-Compton 
(SSC) emission to explain the different epoch detections 
by EGRET and HESS, respectively. 

To explain both, the observed spectrum of LS 5039 
from radio to VHE 7-rays and the variability, we have 
applied a broad-band emission model, based on a freely 
expanding magnetized jet dynamically dominated by cold 
protons and radiatively dominated by relativistic leptons 
(Bosch- Ramon et al. 2005cjl . In this paper we present the 
results obtained after applying such model to LS 5039. 
The compiled observational data covering the full electro- 
magnetic spectrum is given in Sect. [5] In Sect. [3] a brief 
description of the model is provided, and the parameters of 
the model are discussed. The results obtained for LS 5039 
are presented in Sect. 0J and the work is summarized in 
Sect. 

2. Observational data of LS 5039 

We compile here the observational data from the lit- 
erature. The radio data, obtained at 20, 6, 3.5 and 
2 cm wavelengths with the VLA, are from Marti 
et al. I|1998|) . The optical and infrared data, covering the 
bands UBVRIJHKL, are from Clark et al. p1(TT|l and 
Drilling (199lJ, and have been corrected for absorption 
{A v = 4.07, Casares et al. IMiBl . The X-ray data (3- 
30 keV), obtained with RXTE, are from Bosch-Ramon 
et al. <|2005a|) . being diffuse background subtracted. The 
(20-430 keV) data, obtained with the BATSE Earth oc- 
cultation technique, are from Harmon et al. i|2004[) . The 
(1-30 MeV) data are from Zhang et al. dHU), obtained 
with COMPTEL, whose flux values, however, should be 
taken as an upper limit. This is due to the source region 
defined by the COMPTEL 7-ray source, GRO J1823-12, 
at galactic coordinates (1=17. 5°, b = —0.5°), contains 
in addition to LS 5039/3EG J1824-1514 two additional 



EGRET sources (Collmar et al. l2UU31 Zhang et al. !2UU5|l . 
The high-energy 7-ray data have been taken from the third 
EGRET catalog of high-energy 7-ray sources (>100 MeV; 
Hartman et al- E^jl . The VHE 7-ray (MOO GeV) data, 
obtained with HESS, are from Aharonian et al. (2005a). 
The values of the magnitudes and fluxes have been trans- 
formed in luminosities using the updated value for the 
distance of 2.5 kpc (Casares et al. 2005). All these obser- 
vational data are plotted in Fig. 

2.1. X-ray lightcurve 

The X-ray lightcurve of LS 5039 observed by RXTE, 
folded in phase with the new orbital period value, shows a 
clear smooth peak at phase 0.8 and another marginal one 
at phase 0.3 (Bosch-Ramon et al. 2005a). This lightcurve is 
hardly explained by accretion through a spherically sym- 
metric wind of a typical O type star, due to the high ve- 
locities of the wind at infinity, and also because of the con- 
straints on the stellar mass loss rate (^ 10~ 6 M yr , 
Casares et al. I2005[) . It has been proposed that certain 
asymmetries in the wind could produce these peaks in the 
X-ray lightcurve (Bosch-Ramon et al. 2005a). 

3. Modeling 

3.1. An accretion model for LS 5039 

To reproduce approximately the variations in the X-rays, 
we have assumed the existence of a relatively slow equa- 
torial wind (outflowing velocity ~ 6 x 10 7 cm s" 1 ) that 
would increase the accretion rate up to the (moderate) 
levels required by the observed emission in the context 
of our model, and would produce a peak after perias- 
tron. Although the existence of low velocity flows is not 
yet well established, it seems that it would not be a rare 
phenomenon among the massive wind accreting X-ray bi- 
naries. This is, for instance, the case of the 09.5V star 
BD +53°2790, whose ultraviolet spectrum reveals an ab- 
normally slow wind velocity for its spectral type (Ribo 
et al. 12005). Only a small fraction of the wind being 
slower at the equator is enough to increase significantly 
the accretion rate along the orbit if one assumes pseudo- 
synchronization. In this case, the orbital plane and the 
plane perpendicular to companion star rotation axis will 
be likely the same (Casares et al. I2005fl . allowing for the 
compact object to accrete from this slow equatorial wind 
all along the orbit. 

To explain the second peak we have introduced a 
stream, with a velocity of 5 x 10 6 cm s _1 and an additional 
mass loss rate of ~ 1% of the normal wind, that is formed 
near periastron passage when the attraction force between 
the two components of the system is at its maximum. This 
slower matter stream would produce a very delayed peak 
at phase 0.7. The formation of a matter stream, due to 
gravitational stress at the smallest orbital distances, is 
not unusual in close binary systems like LS 5039. We have 
used the approach followed by Leahy l(2002|l to explain 
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Fig. 1. Accretion curve along an orbital period inferred 
from the observed X-ray lightcurve in Bosch-Ramon 
et al. (2005a), adopting the accretion model developed for 
this work. 

the X-ray lightcurve of GX 302-1, also with a peak before 
periastron passage. Although there is still a shift of 0.1 in 
phase between the model peaks and the observed ones, it 
is important to note that there is likely an accretion disk, 
and the matter that reaches the disk spends some time in 
it until part of this matter is ejected as a jet (a lower-limit 
for the accretion time is the free fall time, implying a shift 
of ~0.06 in phase, when falling from 10 12 cm). The ac- 
cretion rate evolution obtained from the model described 
here is shown in Fig. ^ 

3.2. A leptonic jet model for LS 5039 broad-band 
emission 

A new model based on a freely expanding magnetized 
jet, whose internal energy is dominated by a cold proton 
plasma extracted from the accretion disk, has been de- 
veloped in a previous work (Bosch-Ramon et al. 2005c). 
The cold proton plasma and its attached magnetic field, 
that is frozen to matter and below equipartition, provides 
a framework where internal shocks accelerate a fraction of 
the lcptons up to very high energies. These accelerated lep- 
tons radiate by synchrotron, relativistic Bremsstrahlung 
with the stellar wind (external) and within the jet (inter- 
nal), and IC processes. In this model, the external seed 
photons that interact with the leptons of the jet by IC 
scattering come from the star, the disk and the corona. A 
blackbody spectrum is assumed for the star and the disk, 
and a power-law plus an exponential cut-off spectrum for 
the corona. The SSC radiation has also been computed as 
well as the Bremsstrahlung and Compton sclf-Compton 
radiation, being the radiation of the last two mechanisms 
negligible. 

The dissipated shock kinetic energy that goes to rela- 
tivistic particles comes from the mean bulk motion kinetic 



energy, directly related to the kinetic energy carried by 
shocks in the plasma. The number of relativistic particles 
that can be produced along the jet is constrained by the 
limited capability of transferring energy from the shock 
itself to the particles, and by the fact that the relativistic 
pressure must be kept below the cold proton pressure. The 
particle maximum energy is limited by the acceleration ef- 
ficiency and the local energy losses in the available space. 
The acceleration efficiency depends on shock properties: 
shock strength, shock velocity, magnetic field direction, 
and diffusion coefficient (e.g., see Protheroe ll999[l . It has 
been parametrized due to the lack of knowledge of the 
specific details of the shock physics for this case. 

In our scenario, the jet formation is mainly constrained 
by a characteristic launching radius, where the ejected 
matter gets extra kinetic energy from the accretion reser- 
voir, and by the accretion energy reservoir itself, which 
varies along the orbit. The modeling of the accretion rate 
changes allows to introduce variability in a consistent way. 
Moreover, since the interaction angle between jet leptons 
and star photons changes with the orbital phase, variabil- 
ity naturally appears at very high energies (see below). 
Finally, we have approximately computed the effects pro- 
duced on the SED by the photon absorption under the 
external photon fields as well as the IC radiation of the 
first generation of secondaries inside the binary system. 
In Fig. [3 we show the photon absorption coefficient as a 
function of the photon energy, at different distances from 
the compact object, during periastron and apastron pas- 
sages. It is seen that the corona and disk photon fields 
(regions of ~ 10 s cm) attenuate slightly the high-energy 
7-ray photons in the inner part of the jet. At VHE, star 
photon absorption dominates for distances ^ 10 13 cm, i.e., 
within the binary system. At greater distances, absorp- 
tion goes down to negligible values. Photon absorption is 
higher at periastron than at apastron passage, due to the 
dependence of the stellar photon density with the orbital 
distance, although for large distances from the compact 
object, the opacities at both orbital phases become equal. 
Our results are similar to those obtained by Bottcher & 
Dermcr ( 2005) and Dubus H2005fl . that present a detailed 
analysis of the effect of 77 absorption of VHE 7-rays near 
the base of the jet of LS 5039. Moreover, the energy de- 
pendence of the opacity on the ambient photon fields to 
the 7-ray propagation and its variation with the distance 
is not much different from that obtained for the micro- 
quasar LS I +61 303 at its periastron passage (Romero et 

al. 12003 • 

Finally, we note that our model does not take into ac- 
count the radiative contribution from relativistic protons 
that could be accelerated along with the leptons. These 
protons can cool through interactions with ions from the 
stellar wind producing additional 7-ray emission at high 
energies (Romero et al. 120031 12005). 
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Table 1. Primary parameter values. 



Parameter and unit 



value 



e: eccentricity 0.35" 

a: orbital semi-major axis [cm] 2.2 x '. 

8: jet viewing angle [°] 24. 9 a 

m w : stellar mass loss rate [M© yr -1 ] ~7x! 

M x : compact object mass [Mq ] 3.7 a 

i?*: stellar radius [i? Q ] 9.3" 

Mi,: stellar mass [M ] 22.9 a 
L t : stellar bolometric luminosity [erg s _1 ] 7x10 

T*: stellar surface temperature [K] 3.9 x 

p: electron power-law index 2.2 6 

Tjct: jet injection Lorentz factor l.l c 

y: jet semi-opening angle tangent 0.1 c 



log (photon energy [eV]) 



Fig. 2. The photon absorption coefficient as a function of 
the energy for photons emitted in the jet at different dis- 
tances from the compact object, at periastron and apas- 
tron passages. 



t 12 a 
. — 7 a 



38 a 

10 4 a 



kTdisk'- disk inner part temperature [keV] 0.1 

p C o r : corona photon index 1.6 

zq: Jet initial point [i?Sch] 50 

k: jet-accretion rate parameter 0.1' 



a Casares et al. 1200511 . 'Bosch-Ramon et al. 12005aL c Paredes 
et al. (1217021 . d Fender J2T70H . Hujeirat £B3M . 



3.3. Parameters of the model 

The model involves a large number of parameters because 
it aims at reproducing the broad-band spectrum by tak- 
ing into account several physical mechanisms. However, in 
the case of LS 5039 an important fraction of the param- 
eters are well established from observational data or can 
be easily estimated from some reasonable assumptions. We 
call these parameters 'primary parameters'. The remain- 
ing 'secondary' parameters can be reduced to a manage- 
able number. In Table^we present the primary parameter 
values of our model. The first half of them are the param- 
eters describing the binary system and their components. 
We adopt updated values taken from Casares et al. ( 2005 ) . 
The electron power-law index (inferred from X-ray pho- 
ton indices found by Bosch-Ramon et al. I2()()5a|l . the jet 
ejection Lorentz factor (Paredes et al. 12002(1 and the tan- 
gent of the jet semi-opening angle (Paredes et al. I2UU2[) 
are parameters estimated and/or adopted previously by 
other authors. Regarding the jet ejection Lorentz factor, 
although it seems well established that the jet is mildly 
relativistic (Paredes et al. 2002), no direct jet velocity es- 
timate has been obtained until now and we have fixed the 
Lorentz factor to 1.1. The values adopted for the disk inner 
part temperature and the corona photon index correspond 
to typical values found in other microquasars, because the 
X-ray spectrum of LS 5039 only provides clear upper limits 
to the luminosities of both disk (Ldisk — 10 34 erg s _1 ) and 
corona (Lcor ~ 10 34 erg s" 1 ) (Bosch-Ramon et al. 2005a). 
The jet-accretion rate parameter, n, that gives the ratio 
between the jet matter injection rate and the accretion 
rate, has been taken to be similar to observational and the- 
oretical estimates (see, e.g., Fender l2"()01l Hujeirat 2004). 

Next, we describe the secondary parameters of the 
model, which have been adopted in order to better re- 



Table 2. Secondary parameter values. 



Parameter and unit 


value 


£: shock energy dissipation efficiency 


0.5 


/b: B to cold matter energy density ratio 


0.01 


rj: acceleration efficiency 


0.1 


V,,: equatorial wind velocity [cm s" 1 ] 


6 x 10 7 


rhj t: jet injection rate [Mq yr _1 ] 


4 x 10~ 10 


max. ratio hot to cold leptons 


0.1 


r\: launching radius [J?s c h] 


9 


Vstroam: periastron stream velocity [cm s ] 


5 x 10 6 



produce different spectral and variability properties of the 
source. They are summarized in Table The shock en- 
ergy dissipation efficiency, £, that gives the energy dissi- 
pated by shocks in the jet, has been taken (Bosch-Ramon 
et al. I2()05c|l to be at most the half of the kinetic energy 
of the shock to get enough luminosity at energies beyond 
X-rays 1 . The ratio of the magnetic field (B) energy den- 
sity to the cold matter energy density has been fixed to 
Jb = 0.01 taking into account X-ray flux constraints from 
observations (Bosch- Ramon et al. I2005afl . The accelera- 
tion efficiency parameter, r/ (xq c Bc), has been fixed such 
that the maximum energy photons reach the energies at 
which the source has been detected by HESS. This value 

1 This does not mean that half of the shock kinetic energy 
is dissipated through shocks, but that it might be so in case 
that the energy losses would require it, i.e., the shock could 
be radiatively efficient up to a 50% of the kinetic energy in 
the regions where energy losses are very strong (the total jet 
energy lost through radiation is about 15%). 
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for r\ is in agreement with that of a strong and transrela- 
tivistic shock with a magnetic field parallel to its direction 
of motion and a diffusion coefficient close to the Bohm 
limit (Protheroe 1999J. The stellar wind velocity in the 
vicinity of the compact object (K,), estimated in Sect. [3 
is likely affected by peculiarities of the star (e.g., close to 
but not suffering Roche-lobe overflow, fast rotation, etc; 
see Casares et al. I2005|) . This velocity implies, along with 
the known stellar mass loss rate (see Table ^) and assum- 
ing isotropic flux, a mean accretion rate along the orbit 
of m acc ~4x 1(T 9 M Q yr" 1 (i.e. - 5% of the Eddington 
rate). From this value and k, the jet matter injection rate 
can be derived: m ]C t ~4x 10~ 10 Mq yr _1 . In order to try 
to satisfy both, the assumption that the jet is cold matter 
dominated and the observed optically thin nature of the 
radio spectrum, we have fixed the maximum number of 
relativistic leptons 2 to be £ = 0.1 of the total number of 
particles in the jet. This is higher than the value used by 
Bosch- Ramon et al. (2005c). The implications of this are 
discussed in Sect. 0] Finally, the launching radius can be 
estimated to be ~ 9i?sch, since the energy to eject the jet is 
assumed to be accretion energy transfered and the Lorentz 
factor is roughly known (Bosch-Ramon et al. [2005c ) . 

We note that we are not fitting data but reproducing 
a broadband spectrum as a whole. This implies that the 
values of the parameters used by the model are approx- 
imated. The parameters in Tables ^ and d n °t obtained 
directly from observations have a range of validity within 
10% for p, r jct , kT disk , p COI , V w and Stream! an d within a 
factor of a few for zo, K, £, fs, i}, ^jetj C an d r \- 

4. Results 

4.1. Spectral energy distribution of LS 5039 

Figure shows the computed SED (continuous thick line) 
of LS 5039 at the periastron passage (phase 0.0) using 
the physical parameters listed in Tables Q]and[5] We have 
plotted also the contributions from the star, disk, corona 
and synchrotron radiation as well as the inverse Compton 
emission and relativistic Bremsstrahlung. In the same fig- 
ure we plot also the observational data. 

There is a general agreement between the model and 
the observed data, although there are some discrepancies 
that require some explanation. First at all, we note that 
the multiwavelength data were not taken simultaneously, 
and this can be an important source of discrepancy for 
any time-dependent model. 

At the radio band, the computed spectrum has a spec- 
tral index of 0.3, instead of the observed ^0.5 (Marti 
et al. 119981 where F v cx z^ -0 - 5 ). It seems that an addi- 
tional process, not contemplated in our scenario, rises up 
the number of electrons emitting optically thin radio emis- 
sion (e.g., interactions with the wind or the environment). 
Concerning the spatial distribution of radio emission along 

2 The number of relativistic leptons is fixed by £, f, and the 
condition that their relativistic pressure cannot dominate the 
cold particle pressure. 



the jet, our computed fluxes at 5 GHz are about 10, 6 and 
0.1 mJy for the inner part of the jet (< 1 AU), at mid- 
dle distances from the origin (1-1000 AU) and beyond (> 
1000 AU), respectively. This roughly reproduces the ex- 
tended emission of the radio source observed using VLBI 
techniques (Paredes et al. [2002 ) . 

At the X-ray band, most of the emission is synchrotron 
radiation coming from the relativistic electrons of the jet 
and only a marginal contribution comes from the disk 
and almost a null contribution from the corona compo- 
nents (see also, e.g., Markoff et al. I2001fl . Although the 
XMM-Newton data does not show evidences of a thermal 
disk component (Martocchia et al. l2005fl . we have adopted 
in our model a disk with moderate luminosity. The role 
played in the model by this component is negligible, and 
its contribution appears in Fig. [3] as a small bump over 
the power-law spectrum below 1 keV. The computed X- 
ray luminosities and photon indices are very close to those 
detected so far (see next section and also Bosch-Ramon 

et ai.Ennsai- 

Gallo et al. i20Q3|) have proposed that radio emitting 
X-ray binaries, when in the low-hard state, follow a partic- 
ular correlation concerning radio and X-ray luminosities. 
In the context of this correlation, LS 5039 is pretty under- 
luminous at X-rays, presenting a permanent low-hard-like 
X-ray state. This could be explained by the fact that its ra- 
dio emission is optically thin, coming mainly from regions 
outside the binary system, instead of being optically thick 
and completely core dominated (as in the case of mod- 
els such as the one presented by Markoff et al. .20011 see 
also Heinz & Sunvaev l2003[) . As noted above, it might be 
qualitatively explained by the increase in relativistic elec- 
trons predicted by our model, although the issue requires 
further investigation since our radio spectrum is still too 
hard. 

At BATSE and COMPTEL energies, the main con- 
tribution still comes from synchrotron emission of the jet. 
There is a small disagreement between the predicted lumi- 
nosity and the one observed by BATSE. However, we note 
that these data are averaged along many days (Harmon 
et al. 12004(1 . whereas the plotted SED is computed at 
phase 0.0, when the accretion is significant though not 
the highest. At the COMPTEL energy range, as noted 
above, the moderate disagreement between the computed 
and observed data is likely produced by the fact that 
the MeV radiation coming from the COMPTEL source 
GRO J1823-12, to which LS 5039 is associated, might be 
a superposition of the MeV emissions of the three differ- 
ent EGRET sources that are in the field of view (Collmar 
et al. l2U05j l. This means that the data from COMPTEL 
should be actually taken as an upper limit. In such a case, 
our computed values would be consistent with the data. 

At high-energy 7-rays (0.1-10 GeV), the computed lu- 
minosities are similar to those observed by EGRET while 
the computed photon indices are slightly harder. This fact 
could be naturally associated with the low temporal reso- 
lution of EGRET, with the shortest viewing periods span- 
ning for one week and providing time averaged spectra 
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log (photon energy [eV]) 

Fig. 3. Comparison of the spectral energy distribution of LS 5039 computed from the present model, at the periastron 
passage and using the parameters of Tabled with the observed data. We show the different components of the emission 
as well as the sum of all of the them, which is attenuated due to 77 absorption (solid line). The points observed are from 
Marti et al. (1998| (VLA, diamond), Clark et al. (2001) and Drilling (1991) (optical, circle, corrected of absorption), 
Bosch-Ramon et al. ^005a|l (RXTE, square), Harmon et al. l|2lM|l (BATSE, diamond filled), Collmar et al. l|2UU5jl 
(COMPTEL, square filled), Hartman et al. (flt^l (EGRET, circle filled) and Aharonian et al. (|2005a[l (HESS, triangle 
down filled). The arrows in the EGRET and HESS data represent upper limits (3a). 



that could mask harder ones. We note however that, be- 
cause of the possible effect of pair creation in the jet, the 
slopes of the SEDs at high-energy 7-rays, unlike the com- 
puted luminosities, are not very reliable (see Bosch-Ramon 
et al. |2D05c). The dominant component at this band of 
the spectrum is star IC, whereas the synchrotron radiation 
is only important at the lower part of EGRET energies. 
The emission of secondary pairs created within the binary 
system is not significant. 

At VHE 7-rays, star IC emission dominates. The spec- 
trum reaches several TeV, and it is limited by the max- 
imum particle energy, which can reach ~ 10 TeV. The 
contribution of the star IC scattering is favored during 



the periastron passage due to a higher photon density as 
well as a smaller interaction angle between the electrons 
and the photons, which influences the IC scattering. The 
contribution from the the SSC component is minor, al- 
though is more significant at apastron than at periastron. 
The VHE spectrum is strongly affected by the 77 opac- 
ity due to pair creation within the binary system, being 
the absorption stronger at the periastron passage. It ap- 
pears therefore unavoidable an important decrease of the 
emission as well as an additional source of variability at 
those energies. Also, a hardening of the emission, as seen 
in Fig. |3 could take place. Our model underestimates the 
TeV emission, although the difference is less than one or- 
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der of magnitude. Possibly, some high energy processes in- 
volving particle acceleration and emission are taking place 
in the jet at middle scales, where the 100 GeV photon ab- 
sorption is not important. This could also be related to 
the optically thin nature of the radio emission. This outer 
emission would not be as variable as the emission from 
the compact jet. In any case, a deeper treatment at this 
energy band concerning VHE emission and pair creation 
phenomena is required (Khangulyan & Aharonian 2005 ) . 

The results obtained by applying other leptonic mod- 
els (e.g. Bosch-Ramon et al. 2005b:, Dermer & Bottcher 
2005) do not differ substantially from the comprehensive 
model presented here, which can explain the broad-band 
emission from LS 5039. Otherwise, hadronic models can 
reproduce the emission at VHE energies (Romero et al. 
2003) though it is hard to explain how typical EGRET 
fluxes could be achieved (Romero et al. l2005|) . Comparing 
both types of model, it seems that the leptonic ones can 
explain better than the hadronic ones the broadband spec- 
trum at least up to ~ 1 GeV, being both suitable at higher 
energies since the energetic requirements are not as strong 
as those to explain the whole emission from the source. 



4.2. Variability properties of LS 5039 

Once the accretion model has been applied to reproduce 
roughly the X-ray lightcurve, we explore what are the con- 
sequences at the other bands of the spectrum. In Fig. 0] 
we show the flux variations along the orbit at four en- 
ergy bands. At radio and X-ray energies, the predicted 
flux variations are well correlated with accretion, since 
they have similar shapes to the accretion orbital curve 
(see Fig. The predicted amplitude variations are of 
about one order of magnitude in radio and a factor of 
three at X-rays, although we note that, if the radio emit- 
ting region is large enough, say 10 15 cm, the variations 
at these wavelengths will be significantly smoothened by 
light crossing time limitations (as it seems to be the case, 
see Ribo (2003J. The variations of the radio spectral in- 
dices along the orbit are not treated here, due to the model 
limitations at radio. However, we mention that the X- 
ray photon index shows a trend similar to that found by 
Bosch- Ramon et al. H2005af) . i.e. higher fluxes when harder 
spectra. Nevertheless, the anticorrelation is not as strong 
as that observed. Emission at EGRET energies varies al- 
most one order of magnitude, although the shape of the 
lightcurve is different because of the higher rate of IC in- 
teraction during phase 0.0 due to angular effects. At VHE, 
however, the emission is strongly absorbed by the stel- 
lar field, producing a dip around periastron passage. At 
phases 0.6-0.9, there is a long standing peak that is related 
to the stream related peak at phase 0.7 (0.8) that is not 
close enough to the stellar companion for banishing due 
to 77 absorption. A similar peak seems to be present in 
HESS data (Aharonian et al. 2005a) when folded in phase 
with the new ephemeris of Casares et al. H2005II , although 
further observations should be carried out to confirm this. 



5 GHz 




34.10 1 

34.40 3 00 MeV 




0.0 0.2 0.4 0.6 0.8 1.0 

orbital phase 

Fig. 4. The panel shows the flux variations along the orbit 
for the four energy bands considered here. We note that 
emission at very high-energy 7-rays is strongly attenuated 
by the photon absorption in the star photon field. 

At phases 0.2-0.3, a smaller peak is seen, which might be 
detected eventually in future VHE observations. We note 
that, as it is seen in Fig. [21 if TeV emission comes from 
distances typically larger than the semi-major orbital axis 
it would not be significantly attenuated by photon absorp- 
tion. In such a case, variability would be limited by the 
size of the emitting region. 

5. Comments and conclusions 

We have applied to the microquasar LS 5039 a cold mat- 
ter dominated jet model that takes into account accretion 
variability, the jet magnetic field, particle acceleration, mi- 
croscopic energy conservation in the jet, and pair creation 
and absorption due to the external and internal photons, 
as well as the emission from the first generation of secon- 
daries under the external photon fields. 

The computed SED from radio to VHE 7-rays has been 
compared with the available observational data, gathered 
at different bands of the electromagnetic spectrum, show- 
ing in general a good agreement. Although the model is 
limited by its phenomenological nature, we can extract 
some general conclusions that could be useful to under- 
stand the physics behind these objects as well as for going 
deeper in more basic theoretical aspects concerning the jet 
plasma characteristics, particle acceleration and jet forma- 
tion phenomena. In particular, in the context of our model, 
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we can state that the jets are radiatively efficient (~ 15%) 
and sites of particle acceleration up to multi-TeV energies. 
Otherwise, accretion processes are radiatively inefficient. 

The strong absorption of VHE 7-rays by the stellar 
photon field seems to make hard the detection of high 
fluxes of radiation if they are generated in the inner re- 
gions of the jet, closest to the compact object. Possibly, 
since the source is not very faint at these energies, some 
high energy processes involving particle acceleration and 
emission are taking place in the jet at middle scales, where 
the 100 GcV photon absorption is not important. This 
could also be related to the high fluxes and optically thin 
nature of the radio emission. Variability at these energy 
bands could provide information about the location of the 
particle acceleration processes and emitting regions. 

Changes in the accretion rate due to the orbital eccen- 
tricity and a particular wind density profile can explain 
the spectral and variability properties observed at X-rays. 
It hints to a companion stellar wind that is not spheri- 
cally symmetric nor fast in the equator, since otherwise 
the 7-ray radiation could not be powered. 

We conclude that, with reasonable values for the dif- 
ferent parameters, our model seems to be good enough as 
to describe the general features of the spectrum of LS 5039 
from radio to TeV 7-ray energies, the leptonic model being 
a reasonable candidate for explaining the 7-ray emission 
on energetic grounds. 
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